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^ , We investigate the fate of particle production in an expanding universe dominated by a perfect fluid 

. with equation of state p = ap. The rate of particle production, using the Bogolioubov coefficients, are 

determined exactly for any value of a in the case of a flat universe. When the strong energy condition 
^ . is satisfied, the rate of particle production decreases as time goes on, in agreement to the fact that the 

four-dimensional curvature decreases with the expansion; the opposite occurs when the strong energy 
' condition is violated. In the phantomic case, the rate of particle production diverges in a finite time. 

^ ' This may lead to a backreaction effect, leading to the avoidance of the big rip singularity, specially if 
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^ ; 1 Introduction 



• rH , The phenomenon of particle production in an expanding universe is due, essentially, to the fact that in 

^ • curved space-time the vacuum is not unique [HH]. As the universe evolves, and the curvature changes, the 



vacuum state also changes: the initial vacuum state, representing the state with no particle, becomes later 
a multiparticle state. The particle production is directly connected with the curvature of the universe. If 
the universe is spatially flat, the cosmic evolution leads asymptotically to a Minkowski space-time, where 
the phenomenon of particle production does not occur anymore. However, this is true only if the strong 
energy condition is verified: the energy density p and the pressure p must satisfy the relation p + 3p > 0. 
If the strong energy condition is violated the particle production should be zero initially, increasing as 
the universe evolves, in connection with the increase of the four-dimensional curvature with time when 
p + 3p < 0. 
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In order to compute the particle production in a given cosmological model, it is necessary to fix an 
initial vacuum state. If the universe is always decelerating, there is no natural and unique choice for 
this, since all modes are initially outside the Hubble radius feeling strongly the curvature of the space- 
time. However, a natural choice for the initial vacuum state can be done for the case of an accelerating 
universe, since all physical modes are initially well inside the Hubble radius, and they behave as in a 
Minkowski space-time. This is one of the reasons why the primordial inflationary scenario is so successful 
theoretically: it encodes naturally a mechanism for quantum fluctuations, which will be later the seeds 
for the large scale structure of the universe [3]. 

Observations indicate that we live today in an inflationary phase, since the universe is accelerating 
[HIS]. Hence, the mechanism of particle production may play again a very important role today. Suppose 
we define a quantum state for a given field today. The question we want to address concerns the rate of 
particle production. In particular, we would like to know if such particle production is so important that it 
can alters, by back-reaction, the later evolution of the universe. The back-reaction due to quantum effects 
has already studied in the case of the cosmological constant [6] . The back-reaction may be particularly 
relevant for the case the universe is dominated by a phantom fluid since, classically, in such a situation 
the universe will inevitably evolve towards a new singular state, called the big rip. We remember that 
the observational data favors somehow a phantom scenario today [7] . A phantom fluid has many special 
features. Local configurations of a phantom fluid leads to regular black holes fS]. In cosmology, a universe 
dominated by a phantom fluid implies a future singularity that will occur in a future finite proper time 

We will solve the rate of the particle production when the universe is flat and dominated by an 
equation of state p =^ ap, for an arbitrary value of a. A simplified scenario, not considering the different 
phase transitions that occur in the real universe, will be adopted. The calculations will be performed 
for a massless scalar field, for which the corresponding Klein-Gordon equation is solved. Through the 
quantization of the field, the rate of particle may be determined using the technique of Bogolioubov's 
coefficients. The general solution for this problem is expressed in terms of Bessel functions. The rate of 
particle production is determined exactly for any value of a. The connection of this rate with the strong 
energy condition is established. It is shown that for a phantom fluid the rate of particle production 
diverges as the big rip is approached. This may lead to a back-reaction effect changing the effective 
equation of state of the universe, implying perhaps the avoidance of the big rip itself: the particle 
produced does not necessarily obeys the phantom equation of state and may become the dominant 
component of the universe if— 1 > a > — |. Hence, the back- reaction of the particle production may 
be an ineluctable mechanism to the avoidance of the big rip if the pressure is not excessively negative. 
The limiting case, a = — |, has already been found before in a complete different context, that is, in the 
evolution of classical scalar perturbations in phantom cosmological models [TOl H] . 
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This paper is organized as follows. In the next section we review the formalism of particle creation 
in a curved space-time. In section 3 we apply the formalism to a perfect fluid cosmological model, and 
we determine the rate of particle creation for any value of a. In section 4 we present our conclusions. 
Even if many steps of the formalism used and of the calculation performed are quite well known, we will 
present them with some details in order to trace the main features of the final results. 

2 Basic formalism of particle creation 

In the Minkowski spacetime, the Lorentz invariancc allows to identify a unique vacuum state. However, 
in a curved space-time generally the Lorentz symmetries are not present. As a consequence, there is no 
unique vacuum state, and the particle concept (as excitation of a vacuum state) becomes ambiguous. 
The nonuncity of the vacuum state is responsible for the mechanism of " particle production" , which is 
due, in some sense, to the fact that different observers define different states as the "vacuum state". The 
way this mechanism operates is quite model dependent. A general way to analyse the problem is to 
consider a quantum field "living" in a given space-time, and to obtain the particle production via the 
Bogolioubov transformation that connects different vacuum states in a given model. In the present work 
we will consider a free, massless scalar field "living" in a flat Friedmann- Robertson- Walker space-time 
whose dynamics is determined by a barotropic fluid such that p — ap. We will consider a as an arbitrary 
parameter. In particular, a can be negative. 

The Lagrangian density describing a massive scalar field tp non-minimally coupled to gravity is given 

by 

where R is the Ricci scalar and ^ is a coupling constant. The so-called minimal-coupling is obtained when 
^ = 0, while the conformal coupling requires C = ^- We work in natural unities such that G = c = H = 1. 
This Lagrangian density implies the following equation of motion: 

D(p + m^ip + ^R.ip = (2) 

The scalar product of a pair of solutions of the Klein-Gordon equation ([2]) is defined by: 

(^ipi{x),^2{x)^ (^ipUx)X Mx)^ [-.gs(a;)]^dS^ (3) 

where c?E'^ = (En^ with dS being the volume element of the spatial-type hypersurface E, and where 
rt^ is a timelike unitary vector normal to that hypersurface. This scalar product is independent of the 
choice of the hypersurface. 

We proceed by quantizing the scalar field in the canonical way. The conjugate momentum is defined 
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by 

where the dot means derivation with respect to time. The field variable and its conjugate momentum 
satisfy the usual commutation relations: 

[^{t,x),v{t,x')] = , 

[K{t,x),n{t,x')] = , (5) 
[ip{t,x),T:{t,x')] = i5^{x-x') . 

Let us consider {fj} the complete ensemble of positive solutions of equation (l2|); so, {Vj} will be the 
complete ensemble of negative solutions of the same equation. The field may be written under the 
form, 

3 

where Cj and cj are the annhilation and creation operators such that 

[cj,c]]=%.. (7) 

The vacuum state |0) is defined through the condition |0) = 0. In a fiat space-time, the positive 
solutions are given by positive frequencies, and in this way, considering t as the time coordinate, a unique 
vacuum state is obtained, the Minkowski vacuum. 

In a curved space-time, the situation is different. In general the choice of {</?j} is not unique, and as 
consequence the notion of a vacuum state is not unique either. This means that it is not possible to have 
a universal state describing the absence of particle; the notion of particle itself becomes ambiguous. 

Let us consider an asymptotically flat space-time in the past infinity and in the future infinity, but 
not in the intermediate region. The positive frequency solutions in the past infinity are denoted by ipj, 
while the positive frequency solution in the future infinity are denoted by <f)j. These solutions obey the 
following conditions: 

{iPj,Lpj,) = {(j)j,4)j,)^5jj. , (8) 

((p/,^;,) - (^;,<^;,) = -%' , (9) 
(vp,,^;,) = (0„ 0/0 = . (10) 

The field ip may be expressed, at same time, in terms of the functions {fj} or in terms of the function 
{4>j} , since both basis are complete. Hence, 

<^ = E (^^-^^ + 4^3) = E {^^^^ + 4^1) ■ (11) 
3 3 
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The annhilation and creation operators in the past infinity are given by Cj and cj, while bj and ft] are 
the corresponding operators in the future infinity. The vacuum state in the past infinity is defined by 
Cj |0)^^ = Vj, describing the situation where no particle is present initially. The vacuum state in the 
future infinity is defined by hj 10)^^^ = Vj , describing a situation where no particle is present in the 
future infinity. 

The two ensemble of annhilation and creation operators are connected by the relations 

= Y^ioikjbk + Pkjbl) , (12) 

k 

c]=E(/3fei''fe + «fc,4) , (13) 

k 

bk = ^j{aljcj-(3ljc\) , (14) 
bi = '^j{(Xkjc]- PkjCj) . (15) 

The above relations define the Bogolioubov transformation, where ajk and l3jk arc the Bogolioubov's 
coefficients. Since the asymptotic basis are complete, the in-modes can be expressed in terms of the 
out-mode and vice-versa: 

'fij = X] i'^jk^k + f3jk(t>k) . (16) 

k 

'Pk=Yl i^jk'fij - (^ikfj) ■ (17) 
j 

In order the formalism above to be consistent, the Bogolioubov's coefficients must obey the following 
relations: 

= Skr , (18) 

j 

T.i^'hf^n'-Pjka*^') = . (19) 
j 

Now, we can describe the phenomenon of particle creation by a time-dependent gravitational field. 
Using the Heisenberg representation, |0)^^ is the state system at each moment. The number operator 
counting the number of particles in the future infinity is Nk = So, the number of particles in the 

future infinity is given by: 

(A^fc) = (o,in\bibk 0,m^ 
j 

A non null (3jk implies particle creation. 
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3 Particle creation in perfect fluid cosmological models 

The Robertson- Walker metric describing a four dimensional space-time flat, written in terms of the 
conformal time, is: 

ds^ ^a^{T])[dr]'^ ~dx^ -dy^ -dz"^] . (21) 
In terms of the conformal time, the Einstein's equations take the following form: 



pa? , (22) 



{a'Y _ SttG ^ 2 



a 

where p and p are the energy density and the pressure, respectively. We suppose that they are related 
by a barotropic equation of state such that 

p = ap . (24) 

The fact that the sound velocity in this fluid must be equal or smaller than the velocity of the light 
implies that a < 1. From the equations of motion (|22l23p . we obtain 

a{7j)^ao\vf , /?=T-Ar " (25) 
1 + Sa 

If a > —1/3, the times evolves such that < 77 < 00, while, for a < —1/3, the interval is —00 < 
?7 < 0_ This class of cosmological solutions does not interpolate two minkowskian space-time: the four 
dimensional curvature diverges for 77 — > and goes to zero for 77 — > cx3. In the region the curvature goes 
to zero the modes are well inside the horizon and behave as in the Minkowski space-time. This vacuum is 
called Bunch-Davies vacuum and it plays a fundamental role in the origin of quantum fluctuations in the 
inflationary scenarios for structure formation. Fixed the initial vacuum state, the formalism described 
in the preceding section applies for the flat FRW cosmological models and the temporal evolution of the 
occupation number operator can be computed. 

Let us write the solutions of the Klein-Gordon equation under the form 



where the functions ipk (j], x), through a separation of variable, have the form 

^,{rj,x) = {2nyK^''-'^ (27) 
Using the flat Robertson- Walker metric, the massless Klein-Gordon equation takes the form, 

^'^ + 2-ipk + k\k = . (28) 
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Inserting ([27]) into ([28]), we obtain 



a 



Xk"{v)+{k~ — ]xk{v)=0 • (29) 



Using the expression for the scale factor, the above equation reduces to 

Xfe"(^/)+(fc'-^^^)xfc(^)=0 . (30) 

This is essentiahy the same equation governing the evolution of gravitational waves in an expanding 
universe [T^ . 

This equation admits solutions under the form of Hankel's functions: 



Xkiv) = V' 



A,ijW(fc|ry|) + Mf (fchi) 



9=1^-/31 , (31) 



('121 

Hj^ ' ' are the Hankel's functions, {Ak, B^) are constants whose values are fixed by normalizing the 
corresponding modes. To do this, we use the Wronskian relation 

zH(^^ (z)a,i? (z) - zi? (z)9.ilf ) (z) = - . (32) 

Imposing the orthonormalisation of the modes one obtains: 

\A,\'^\B,\'^^ . (33) 

The function Xk (fl) becomes (up to an arbitrary phase) 

Xk{v)^^(e,H('Kk\ri\) + e2H('Hk\v\?j , £1,2 = 0, 1 . (34) 

This solution has been already studied previously (see |13|). but in view of a primordial cosmological 
models, in special in de Sitter or quasi-de Sitter phase. 

With the aim of determining the particle production for the cosmological model, we will rewrite the 
Bogolioubov's transformations in a convenient way. For a minimal coupling and a massless field, the 
Lagrangian density takes the form: 

/:(x) = i[-g(a:)]^{.g^Xa;)^^^(a;)_J . (35) 

With the background described above, the Lagrangian can be rewritten as 

The conjugate momentum is 
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Hence, the Hamiltonian is 



n = l\p^ + + 2-xp 

2 a 



After quantization, the operators fi and p satisfy the commutation relation 

[x,P]=i ■ 



(38) 



(39) 



We can now determine the expressions for the creation and annhilation operators, and b respectively, 
with respect to the Hamiltonian (|38p . To do this, we fix: 



1 



X ■ 



(b + b^) 



P - 



~i\l^{b-b^) 



These operators satisfy the commutation relation 
The Hamiltonian takes the form 



H = kb^b + G{r^)b^ + a*(ry)62 + (fc + 'l.) 

2 a 



where the coupling function is given by cr(?]) = ia' /2a. Using the normal ordering 

■.n:=n- {Q\'H\Q) , 

one obtains 

H = A:6^6 + cr(77)6t%cr*(7y)62 , 



(40) 



(41) 



(42) 



(43) 



(44) 



We will use the Heisenberg representation for which the operators evolve with time and the quantum 
states remain fix. Since there are two independent modes, the Hamiltonian (|44p can be written as a sum 
of two Hamiltonians Tii and 7^2 with the same frequency: 



Ul = kb\bl + (7(7?)6f + CJ*{T^)b\ 

U2 = kb\b2 + cT{jf)bf + cT*{ri)bl 



Thus, the Hamiltonian becomes 



n = k 



b\bi + blb^ + airf) 



We now define two annhilation operators with momenta k and —k: 

bi — ib2 bi + ib2 

The Hamiltonian may be now rewritten in terms of these new operators: 

n = kcicf: + kcl-c^f: + 2a{7j)cicl. + 2a*{f])c^c 



(45) 
(46) 

(47) 

(48) 

(49) 



The new operators obey the foUowing equations: 



del 



drj 



ci.H 



Usmg ([49]), we find: 



del 



whose solutions can be written as 



d?7 



a *■ 



(50) 
(51) 

(52) 



(53) 
(54) 



where c^(0) and ct(0) are the initial values of the operators Cj:{ri) and ct(?7), respectively. The function 
Vk is nothing else than the Bogolioubov coefficient /3ij defined before. 
From (|53|) . evaluate at 77 = 0, one obtains: 



c^{0)=Uk{0)c^{0)+Vk{0)cl.{0) , 



resulting that Mfe(O) = 1 and Wfe(O) = 0. Due to the fact that ci:{ri), ct[ri) 



1, one can write 



resulting that 



\uk{v)f - \vk{v)\^ = 1 



That is the normalization condition. We find also the relations 

duk{r]) 



drj 
dvk{vi) 



-ikukiv) + —vlirj) , 



= -ikvkirj) + —ulirj) 
dij a 



The sum {uk + Vk) satisfies the relation 

(uk + vl) 

Moreover, 



(uk + vl) = 



leading to 



Xk{v) = Ukiv)+v*kiv) 
13 



1 

2 ^*2fcr; 



Xkiv) - ^Xkiv) 



(55) 

(56) 
(57) 

(58) 
(59) 

(60) 

(61) 
(62) 
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Using the solutions (p4| for the field Xfc(f]), with £2 = 0, we find 

vtM = ^^e^'" k^Hfchl) - ^Hj^'Mm)] ■ (63) 



Hence, 



^ (H^,'HkM)H^'\ik\v\) - i^;^\(fc|r;|)fff (fchi) 



(64) 



This expression gives the rate of particle production. It has the following asymptotic behaviours: 



TT_l'k\r]\ 



3q-2 



, 9<i (a>i) ^ . (65) 

fc|ry|^oo , Vg ^ |i;fc(ry)|2 ^ . (68) 

For a > —1/3, 77 ^ implies t ^ 0, while for a < —1/3, 77 ^ implies t 00. The asymptotic 
behavior above just show that, for a > —1/3, there is no natural initial vacuum state, except for the 
radiative case a — 1/3: the initial state is fixed arbitrarily by hand. That is, an decelerating universe 
requires a previous accelerated phase in order to establishes a natural vacuum initial state. For a < —1/3 
there is, in opposition, a natural initial vacuum state, and the particle production is initially zero, 
becoming infinite asymptotically. 

3(l + a) 

In general, the Hubble radius is given by duiri) oc 77 1+3° . The Hubble radius increase with time for 
a > — 1, remaining constant for a = — 1, and it decreases for a < — 1. Hence, in the phantom case, all 
modes are well inside the Hubble radius if we go back far enough in the past, and it is quite reasonable 
to assume that the initial state is the vacuum state Q. In figure 1 the behavior of the scale factor, of 
the Hubble radius and of an arbitrary physical model are displayed for a phantom cosmological model, 
illustrating the situation we have just described. 

For the phantom scenario the future asymptotic occurs in a finite proper time. Hence, the particle 
production may change the final state; in particular, it may lead to the avoidance of the big rip. This 
may happens since the equation of state of the particles created, for a massless scalar field, is given by 
Ps = Ps, and we can expect that the phantom fluid will not become the dominant component anymore. 
In fact, the energy density of perfect fluid scales as pph oc a^^^^+"^ oc rj" 1+3° . If we suppose the energy 



*In the real universe we should perhaps prefer to consider an initial non vacuum state, since there was already a particle 
production in the previous phases of the universe. But, an initial non vacuum state does not change our main conclusions, 
except for some very special cases. 
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Figure 1: Behavior, as function of time, of the scale factor (red), of the Hubble radius (green) and of a typical physical 
mode (violet) in a phantom model. 



density of the scalar particle created scales according to the rate evaluated above, we have ps oc Tyi+a^. 
So the ratio pph/ Ps goes to zero if —1 > a> —5/3 and to infinity if a < —5/3: the back-reaction can be 
effective if the equation of state of the phantom fluid is not deeply negative. 

The particular case a = — must be studied separately. In terms of the conformal time, the solution 
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for this case takes the form 

a{r]) = aoe^'' . 

where ao and A are constants. The Klein-Gordon equation now reads, 

Xk"{v)+D\k{v) = , 



with 



(69) 



(70) 



(71) 



There is no propagation of the quantum modes if D'^ < 0. However, for £) ^ > 0, the solution becomes 



The normalisation implies that 



Using the last expression, we obtain 



Xkiv) = Xoe 



Xo = {2D)- 



-iDr} 



*( \ 11 



, D .A 



(72) 



(73) 



(74) 



Particles, for this case, arc created with a constant rate during all the evolution of the universe for the 
modes that arc inside the Hubble radius, while no particle is created for the modes that arc outside the 
Hubble radius. Note that, in this case, a mode that is initially inside the Hubble radius remains always 
inside it. 



4 Conclusions 

We have evaluate the rate of creation of massless scalar particle in a universe dominated by a perfect 
fluid whose equation of state is given by p = ap. An analytic expression has been found in terms of 
Hankel's functions. Since we have not considered a complete cosmological scenario, with a sequence of 
difi'erent phases, with an initial inflationary phase (as in standard cosmological model), the calculation 
performed here makes sense, strictly speaking, only when the strong energy condition is violated, that is, 
a < —1/3. For such a case, there is a natural initial vacuum state, from which the particle occupation 
number can be determined. However, we formally extended the calculation for any value of a. It has 
been found that for those "dark energy" scenarios the rate of particle creation diverges as f ^ oo, that 
is, in the future infinity. 

We have payed special attention to the phantom scenario. The reason is that a universe dominated 
by a phantom fluid may develop a singularity in a finite future time, the so-called big rip. In this case, 
we have found that it is possible that the energy density associated to the particles created (which obey 
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in present case an equation of state of the type Ps — Ps) becomes dominant over the phantom fluid if 
— 1 > a > —5/3. Hence, the big rip can be avoided if the pressure is not deeply negative. 

It is interesting to remark that a similar critical point, a = —5/3, has been found in the case of 
classical scalar perturbation [ll . In that case, however, the evolution of scalar perturbation may destroy 
the conditions of homogeneity (necessary for the big rip) if a < —5/3, that is, if the pressure is negative 
enough. The result found here is exactly the opposite: quantum effects can be operative in the sense 
of destroying the conditions for the big rip if the pressure is not negative enough. It must be stressed, 
however, that the evaluation made in the present work must be complemented by a study of more general 
quantum fields and by a deeper thermodynamical analysis of the energy balance between the phantom 
fluid and the created particles as the big rip is approached. 
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